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Abstract: We calculate the cross section for ccg-pro duct ion in diffractive DIS with 
finite quark masses at zero momentum transfer t. The calculation is done in the 
leading log(l/xp) approximation and is valid in the region of high diffractive masses 
Q\ ', M (small (3). We apply our cross section formula including both cc- and ccg in a 

Monte Carlo simulation to diffractive meson production at HERA. We compare 
1 our predictions to results of HI using three parameterizations for the unintegrated 

gluon density. 

O 

Ph! 1 Introduction 

In the process of diffractive deep inelastic scattering, 7* + p — > p + X, one can separate per- 
turbative and non-perturbative contributions by filtering out particular diffractive final states. 
Examples of diffractive states which are perturbatively calculable are longitudinal vector parti- 
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cles or final states which consist of hard jets (and no soft remnant). In the latter case the hard 
scale which allows the use of pQCD is provided by the large transverse momenta of the jets, 
and the Pomeron exchange is modeled by the unintegrated gluon density. Another particularly 
interesting example is diffractive charm production, since the charm quark mass justifies pQCD, 
even for not so large transverse momenta of the outgoing quarks and gluons. Calculations for 
the diffractive production of massless open qq states and of massless qqg states have been re- 
ported in [ [I], ^ [!| and in [ ||, respectively, and a comparison of diffractive two-jet and three-jet 
events observed at HERA with these calculations has been presented in [ [| . Final states with 
finite quark masses have been calculated, so far, only for qq production [ [| which is expected 
to be the dominant final state in the region of small diffractive masses (large (3). However, as 
there are recent measurements of diffractive - production from the HI [ |7| and ZEUS [ 
collaborations at HERA, which extend into the small-/3-region, gluon radiation can certainly 
not be neglected, and a full perturbative calculation of qqg is needed. 

In this article we report on a calculation of massive qqg-pro duct ion in diffractive deep inelas- 
tic scattering, and we present a comparison of our cross section formula with the measurements 
of HI [fa. 
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2 Calculation of massive qqg production 



We will follow the study of massless qqg-pro duct ion presented in [|3J]. In particular, we again 
work in the leading-log M 2 approximation, which limits the applicability of our results to the 
small /3-region. Fig. [| shows the notations of the process. As in [ |4j] we restrict ourselves to 
zero momentum transfer, t = r 2 = 0. As usually, Q 2 denotes the virtuality of the photon, 
VW 2 the energy of the photon proton system, M the mass of the diffractive system, and 
x = Q 2 /{Q 2 + W 2 ), y = pq/pl are the Bjorken scaling variables (with I being the momentum 
of the incoming electron). The variable (5 is defined as (3 = Q 2 / (Q 2 + M 2 ), and it is convenient 
to introduce the momentum fraction of the Pomeron by xjp = (Q 2 + M 2 )/(Q 2 + W 2 ). 

We restrict our calculation to the region (leading- log M 2 approximation): 

Q 2 < M 2 < W 2 . (1) 
We use Sudakov variables ki = a^q' + (3(p + k it (with q' = q + xp, kf t = — kf), and we express 
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Figure 1: Kinematics of diffractive qqg production 



the phase space in terms of y, Q 2 , M 2 , m 2 , t, kf, k 2 , with m 2 = m 2 qq + k| {m qq denotes the 
invariant mass of the gg-subsystem) . We obtain the following result: 
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The differential cross sections of 7*p — > ccg + p are given by: 
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and the function T denotes the unintegrated (forward) gluon density^ which is connected with 
the usual gluon density g(x,Q 2 ) through: 



rQ 2 

/ rfl 2 ^(x,l 2 ) ~xg(x,Q 2 ). 
Jo 



(12) 



The ~ sign in the above equation indicates that the relation is valid for large Q 2 . Strictly 
speaking, the kinematics of diffractive qqg production requires the nonforward (skewed) gluon 
density. However, our cross section formula has been derived in the leading- In W 2 , leading- 
lnM 2 approximation, and the use of the gluon density in (12) is valid only in the double 
logarithmic approximation where skewedness is negligeable. 

The parameter ol\ is determined by the on-shell conditions for the final state particles: 



° 1= 2 
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(13) 



1 Note that different definitions for T exist, here we use T as defined in eq.fl). 
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and it varies between and 1. The values of the momenta k 1; k 2 and of m 2 determine the sign 
in eq . (|13"D . The validity of our cross section formula is restricted to the kinematic region where 
the gluon transverse momentum k 2 is not small. 

The quark mass m q enters the calculations in two places. First, the phase space of the 
diffractive system (and so the parameter a\ and the function S) depend upon the quark mass 
via the on-shell conditions for the outgoing particles. Secondly, the propagators of the internal 
fermion lines are modified by a nonzero quark mass which leads to changes in the matrix- 
elements. Apart from the function D(k), eq . (|ll|) , which enters in all four 7*p - cross sections, 
an additional term containing the quark mass emerges in rfo"Jy + (eq.(|3|)). 



3 Comparison with measurements 

Compared to other charmed particles, D* ± mesons are easy to reconstruct which makes them 
attractive objects for testing diffractive charm production. D*^ mesons are identified via the 
decay channel 

D* + -> D\: iow -> (K-n + )n+ ow (and c.c), 

which has a branching ratio of 2.63% [ §]. In the following comparison we concentrate on a 
comparison with the measurement of the HI collaboration [JJ§, who has analyzed data collected 
throughout the years 1995-1997. The amount of data is still quite poor due to the small 
branching ratio of the D* meson, but higher statistics will come from new data. 

We have implemented the cross section formulae for diffractive massive cc production [ || and 
from our expression eq.(Q-|6|) for the massive qqg production into the Monte Carlo generator 
Rapgap [ [10], |n| , which includes full hadronization according to the Lund string model as 
implemented in Jetset/Pythia [ [T2|, |T3|. We have used a fixed strong coupling constant 
a s = 0.25, and a charm quark mass of m q = 1.5 GeV. The transition of the charm quark 
to the D* meson is performed via the Lund heavy quark fragmentation function. Since the 
experimental measurement does not separate the charge of the D* meson, we have added the 
cross sections of both charged D* mesons. 

Events generated with Rapgap are selected within the same kinematic region as in the 
measurement of HI [ [7|], using electron and proton momenta (in the HERA system) of 27.6 
GeV and 820 GeV, respectively, and: 

0.05 <y <0.7 (14) 
2 < Q 2 < 100 GeV 2 
x P < 0.04 
\t\ < 1 GeV 2 

As the measured cross section is integrated over t for \t\ < 1 GeV 2 , we have multiplied our cross 
section formulas (eq.(^-§|)), valid for \t\ = 0, with a phenomeno logically motivated t distribution 
of the form: 
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The D* mesons are experimentally detected in the range: 

IviD^l < 1.5 (16) 
Vt{D*±) > 2GeV 

with the pseudo-rapidity 1] = — logtan(#/2) and the transverse momentum of the D* meson 
measured in the ep laboratory system. 

As can be seen from eqs.(P,|T0|), there is a potential divergency, if the transverse momentum 
of the final state gluon k| approaches zero. In order to avoid the non-perturbative region, we 
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Figure 2: The cross section da/d log 10 xp for diffractive D* production within the kinematic 
range specified in the text. The point are the measured cross section from HI [ [?[. The 
prediction obtained with k cut = 1 (1.5) GeV is shown in a (b). The dashed (dotted) line shows 
the cc (ccg) contribution alone and the solid line is the sum of both. The dGRV unintegrated 
gluon density is used, with ¥ min = 0.4GeV 2 . 



impose a lower cutoff k 2 cut on the gluon transverse momentum. In our calculations we have 
considered k 2c ut = 1 GeV and k 2cn f = 1.5 GeV. For the unintegrated gluon density, which 
enters in eq.(||), we have used three different approaches: the derivative of the NLO GRV [ |I4[ 
gluon density dGRV, the unintegrated gluon density J-(x, l 2 ) obtained in the saturation model 
of Golec-Biernat and Wusthoff [|EJ GBW } and the CCFM [ [T§, |T7|, |I§ unintegrated gluon 
density A(x, l 2 , q 2 ) JS of [ |2lfl , where q defines the evolution scale, related to the maximum 



allowed angle of any emission in the angular ordering approach. 

The unintegrated gluon density can be obtained from the integrated gluon density, if in 
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eq.(|T^) the ~ sign is replaced by an equality sign: 

d xg(x, n 2 



d fi 2 



(17) 

=i 2 



Here we use for xg(x : \i ) the NLO GRV [ [14j] gluon density, since it is the only integrated 



gluon density available, starting at a low value of = 0.4 GeV 2 . Due to the finite Qo in 
any of the available integrated gluon densities, a lower integration limit l^ in ~ Ql in eq.(||) is 
introduced. Variation of this parameter mainly affects the normalization of the cross sections. 
For example, when l^ in is decreased from 1 to 0.5 GeV 2 , the xp-distribution at k 2cu t = 1 GeV 
roughly doubles in the whole xp range. We have chosen to set ^ min as small as it is compatible 
with the definition of the integrated gluon density, so ¥ min = 0.4 GeV 2 for dGRV. The other 
two unintegrated gluon densities, GBW and JS, are defined also for the very small l 2 region, 
and therefore no cut needs to be applied there. In the numerical treatment of the cc production 
cross section neither k 2 ,^ nor l^ in are needed. 

In Fig. H we show the effect of changing the cut k 2 ut on the differential cross section 
da/dlogiQXp, for diffractive D* meson production in the kinematic region specified above and 
compare our prediction to the measurement of HI [ [7]. Also shown in Fig. [| are the indi- 
vidual contributions of cc (dashed histogram) and ccg (dotted histogram). One clearly sees 
the reduction of the ccg contribution when the cutoff in h cut is increased. Since our calcula- 
tion is valid for small xp, we focus only on agreement in the low xp region. A cutoff value 
of k 2 ut = 1.5 GeV 2 seems to be a reasonable choice which we will keep independently of the 
choice of the unintegrated gluon density. With this cut the computed cross section in the lower 
xp bin agrees well with the data. In the upper xp bin, however, the theoretical curve is by a 
factor of about 10 smaller than the data point (ignoring the large error on the measurement). 
In this xp-region it is expected that, apart from Pomeron exchange (which, in our model, is the 
2-gluon exchange) secondary exchanges have to be included: in a perturbative description such 
an exchange corresponds to gg-exchange. Since such a contribution has not yet been included 
into our calculation, it is not surprising that the two-gluon model undershoots the data. 

In Fig. we compare our prediction of the cross section for diffractive D* production as 
a function of log 10 j3 with the measurement of HI (it is understood, that all other variables, in 
particular xp, are integrated over). We also show the individual contributions of cc (dashed 
histogram) and ccg (dotted histogram). Clearly, the ccg is badly needed in order to get closer 
to the data than with cc alone. The measured cross section is slowly rising with decreasing 
log 10 (3. The theoretical curve does not quite follow this rise, the reason for this is the correlation 
between (3 and xp: small (3 values require large diffractive masses M, which due to the kinematic 
restrictions in the analysis, are predominantly produced at large xp. As argued above the large 
xp-region needs secondary exchange which in our approach is not yet included. For illustration, 
we show in Fig. [|& the cross section for diffractive D* production da /dj3. We observe that the 
strong drop in cross section at small (3, as seen in Fig. |3]a is a consequence of plotting da /d log 10 (3 
instead of da /df3. From Fig. |3]& we see that the theoretical curve increases towards low /3, and 
only at (3 < 0.01 is decreasing, which again is a consequence of kinematic correlations. If the 
double differential cross section d 2 a/df3dxp is considered, our prediction shows the expected rise 
towards small (3 at fixed xp. The shape of the theoretical curves as a function of (3 is almost 
independent of the choice of the cutoffs l 2 ^ and k cut , as they mainly affect the overall magnitude 
of the cross section. In the small (3 region, corresponding to large diffractive masses M, also 
the radiation of more than one gluon, such as ccgg need to be considered. Since experimentally 
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Figure 3: The differential cross section da/d\og w (3 for diffractive D* production compared to 
the measured cross section from HI (a). In (b) the differential cross section da/df3 is shown. In 
all cases we use k cut = 1.5 GeV and the dCRF unintegrated gluon density. The dashed (dotted) 
line shows the cc (ccg) contribution alone and the solid line is the sum of both. 



the cross section is defined as e + p — > e' + (D* + X) + Pdiff, where in the diffractive system 
M = D* + X, the hadronic state X is not further specified or measured, multiple soft gluon 
contributions might be present in the data, which have yet not been estimated consistently in 
the perturbative calculations. 

The saturation model of K.Golec-Biernat and M.Wusthoff [ [15|] describes a completely 
different approach to estimate the unintegrated gluon density T{x, l 2 ), which is needed in 
eq.(|]). In this model the total 7*p cross section is described by the interaction of a qq pair 
(dipole) with the proton, and a particular ansatz is made for the dipole cross section. The 
function ^(x, l 2 ) has the form 

Hx ,P ) = gLB> {x) ? e -W* , ^ = ^(^) A/2 , (18) 

and the three parameters of the model are determined by fitting inclusive DIS data (including 
charm with: a = 29.12 mb, A = 0.277, x = 0.41 10" 4 and a s = 0.2 [ 0). For large l 2 , 
F{x, l 2 ) has the meaning of the unintegrated gluon density, but for smaller l 2 the function 
looses this interpretation and has to be viewed as a (model dependent) extrapolation. The 
ansatz eq.(|i~8|) holds for the qq color dipole cross section. The qqg system consists of a color 
triplet, and anti-triplet and a color octet, and one expects that the dominant configuration is 
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Figure 4: The differential cross section da/d log 10 (3 for diffractive production compared to 
the measured cross section from HI (a). In (b) the differential cross section da/df3 is shown. In 
all cases we use k cut = 1.5 GeV and the GB W unintegrated gluon density from the saturation 
model. The dashed (dotted) line shows the cc {ccg) contribution alone and the solid line is the 
sum of both. 



a dipole consisting of two octets: in eq. ([l8|) we therefore rescale the color charge and use: 



Insertion of the ansatz (eqs.( p!8| , |l9D ) into eq.(12) leads to an integrand that vanishes as l 2 goes 
to zero. Therefore, within this model we no longer need any lower cutoff in the l 2 integral, and 
our calculation provides absolute predictions of the cross sections (note, however, that we still 
have the cutoff k cn4 = 1.5 GeV on the final state gluon). Fig. ^ shows the differential cross 
section da/d log 10 xp calculated using eq. ( |T8|JT9| ) compared to the measurement of HI. The 
calculated cross section (cc+ccg) is similar to that of Fig. |3|, and the same discussion applies. 
Also the P distribution (Fig. ^6) is very similar to the previous model. Note that without the 
correcting color factor in the dipole cross section formula the ccg cross section would be larger 
by about a factor of about 1.5. 

The unintegrated gluon density, based on the consistent treatment of color coherence effects, 



is described by the CCFM evolution equation [ 16, 17, 18, 19]. According to the CCFM 
equation, the emission of partons during the initial state cascade is allowed only in an angular- 
ordered region of phase space. In the large (small) x limit, the CCFM equation is equivalent to 
the DGLAP (BFKL) evolution equations, respectively. A solution of the CCFM equation has 
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Figure 5: The differential cross section da/d log 10 (3 for diffractive D* production compared to 
the measured cross section from HI (a). In (b) the differential cross section da/df3 is shown. 
In all cases we use k cut = 1.5 GeV and the CCFM based JS unintegrated gluon density. The 
dashed (dotted) line shows the cc (ccg) contribution alone and the solid line is the sum of both. 



been found, which successfully can be used to describe a bulk of measurements at HERA and 
the Tevatron [|2(| [21], |22|| . However, due to the angular ordering requirement, the unintegrated 
gluon density A(x, l 2 , q 2 ) is now also a function of the evolution scale q, which is related to the 
maximum allowed angle. Here, this scale is set either by the qq pair, or by the final state gluon 
for qqg: 

q 2 = { k f ^ y ! 20 
[ if; for qqg 

with Q t being the vectorial sum of the transverse momenta of the qq pair, and z = (Q 2 + 
m 2 q )/(Q 2 + M 2 ). Since the explicit parameterization of A(x,\ 2 ,q 2 ) from [ |20]| is valid also in 
the very small l 2 region, no cut on l 2 needs to be applied for the integral in eq. The results 
for the differential cross sections as a function of xp and the log 10 (3 shown in Fig. |5]a and b are 
quite similar to those of the unintegrated gluon density from saturation model GBW and or 
from the derivative of the integrated gluon density dGRV. The main difference is the ratio of 
cc and ccg contribution. The enhancement of the (3 distribution by ccg is much stronger than 
in case of the two other models. 
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4 Conclusion 



In this article we have analyzed DIS diffractive charm production (production of D*^ mesons) 
within the perturbative two-gluon model. For the two-gluon amplitude we have used three 
different models: the unintegrated gluon density derived from the integrated DGLAP gluon 
density dGRV, the saturation model of Golec-Biernat and Wiisthoff GBW, and a CCFM-based 
unintegrated gluon density JS (the last two models are parameter-free, the first one depends 
upon a cutoff on the internal momentum integral). In all three cases the calculated cross sections 
are of the same order of magnitude as the data, and within the kinematical region where the 
models apply, the shapes of the cross sections are consistent with the data. Compared to an 
earlier attempt (with the dGRV gluon density) where only cc production had been included in 
the theoretical analysis the present analysis contains, as the new ingredient, also (massive) ccg 
production and leads to a considerable improvement in the agreement with experimental data. 

It is encouraging to see that, for the dGRV gluon density, our analysis of diffractive charm 
production uses the same parameters as in the successful analysis of diffractive jet production 
and we were able to consistently describe both types of processes. 

We view the use of the two-gluon model as part of a more general strategy of analyzing 
DIS diffraction data at HERA. In a first step one would analyze those diffractive final states 
which are dominated by short distances (diffractive jets or states consisting of heavy quarks): 
in these processes the application of the two-gluon model can be justified. In a second step, 
one would try to extrapolate also into kinematic regions where soft physics becomes important. 
In the present analysis, such extrapolations are contained in the saturation model and in the 
CCFM amplitude; the use of the GRV gluon density, on the other hand, requires a momentum 
cutoff. In a final (and future) step one would need to find a QCD-based 'derivation' of the 
extrapolation from hard to soft physics. 

Despite this encouraging success, several improvements in the theoretical part of our 
model should be made. First, since the cross section formula for ccg production has only been 
calculated in the leading log-M 2 approximation, an improvement which extends the 
applicability down to small-M 2 values would be very desirable. For consistency reasons, one 
then would need a NLO-calculation of cc production. Results of such a calculation would also 
allow to eliminate the cut on the transverse momentum of the final state gluon. Next, since 
the region of xp > 0.02 seems to require secondary exchanges, they should be modeled, in the 
framework of perturbative QCD, by qq exchange. Finally, our comparison with data indicates 
the need of ccgg final states: such an extension (at least in the leading log-M 2 approximation) 
should be fairly straightforward. A successful test of the two-gluon model in DIS Diffraction, 
apart form providing a description of charm or jet production at HERA, is also of general 
theoretical interest: the cross section formula for diffractive qq + ng production contains the 
perturbative triple Pomeron vertex which is expected to play a vital role in the unitarization 
of the BFKL approximation. It has been calculated both analytically and numerically, and 
these calculations can be tested experimentally in DIS diffraction dissociation. 

References 

[1] J. Bartels, H. Lotter, M. Wiisthoff, Phys. Lett. B 379 (1996) 239, ERRATUM-ibid. B 
382 (1996) 449, |hep-ph/ 9602363 . 



10 



[2] N. Nikolaev, B. Zakharov, Z. Phys. C 53 (1992) 331. 

[3] E.Gotsman, E.Levin, U.Maor, Nucl. Phys. B 493(1997)354. 

[4] J. Bartels, H. Jung, M. Wusthoff, Eur. Phys. J. C 11 (1999) 111, |hep-ph/ 9903265 

[5] C. Adloff et al, Eur. Phys. J. C 20 (2001) 29, |hep-ex/001205T| . 

[6] H. Lotter, Phys. Lett. B 406 (1997) 171, |hep-ph/96l2415 . 



[7] C. Adloff et al, Phys. Lett. B 520 (2001) 191, |hep-ex/0108047[ 

[8] ZEUS Collaboration, Diffractive D* ± (2010) production in deep inelastic scattering at 
HERA, in Proc. of the 30th International Conference on High Energy ICHEP 00 (Osaka, 
Japan, July, 2000), contributed paper 874. 

[9] Partcile Data Group; C. Caso et al., Eur. Phys. J. C 3(1998) 1. 

[10] H. Jung, Comp. Phys. Comm. 86(1995) 147. 

[11] H. Jung, The RAPGAP Monte Carlo for Deep Inelastic Scattering, version 2.08, Lund 
University, 2002, http://www.quark.lu.se/~hannes/rapgap/. 

[12] T. Sjostrand, Comp. Phys. Comm. 82(1994)74. 

[13] T. Sjostrand et al., Comput. Phys. Commun. 135 (2001) 238, |hep-ph/0010017 . 

[14] M. Gliick, E. Reya, A. Vogt, Z. Phys. C 67(1995)433. 

[15] K. Golec-Biernat, M. Wusthoff, Phys. Rev. D 60(1999) 114023. 

[16] M. Ciafaloni, Nucl. Phys. B 296 (1988) 49. 

[17] S. Catani, F. Fiorani, G. Marchesini, Phys. Lett. B 234(1990)339. 
[18] S. Catani, F. Fiorani, G. Marchesini, Nucl. Phys. B 336(1990)18. 
[19] G. Marchesini, Nucl. Phys. B 445(1995)49. 

[20] H. Jung, G. Salam, Eur. Phys. J. C 19(2001)351, hep-ph/0012143. 

[21] H. Jung, Phys. Rev. D 65(2002)034015, DESY-01-136, |hep-ph/0110034] . 

[22] H. Jung, Unintegrated parton densities applied to heavy quark production in the CCFM 
approach, in Proceedings of the Rinberg workshop on "New trends in HERA physics", 
Ringberg Castle, Tegernsee, Germany. (2001), hep-ph/0109146. 



11 



